Exponential Temperature Dependence of Penetration Depth in Single Crystal MgB2 
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The temperature dependence of the London penetration depth, A(r), was measured in both single 
crystal and polycrystalline MgB2 samples by a high-resolution, radio frequency technique. A clear 
exponential temperature dependence of A(r) was observed at low temperature, indicating s-wave 
pairing. A BCS fit to the lowest temperature data gives an in-plane energy gap A of 30 ± 2 K 
{2/S./Tc — 1.5 ± 0.1), which is significantly smaller than the standard BCS weak coupling value of 
3.5. We find that the data are best described by a two-gap model. 



PACS numbers: 74.25.Nf 

The recent discovery Q of superconductivity at 39 K 
in the binary compound MgB2 has sparked a large num- 
ber of investigations into its physical properties. A cru- 
cial question is whether its high Tc can be explained by 
a phonon mediated pairing interaction within the usual 
BCS-Eliashberg framework. A first step in answering this 
question is to determine the symmetry of the supercon- 
ducting order parameter and the nature of the low energy 
excitations. 

Magnetic penetration depth (A) measurements are a 
powerful probe of the low energy excitations. As A(0) > 
1000 A in MgB2, penetration depth measurements probe 
relatively large distances and are far less sensitive to sur- 
face quality than, for instance, tunneling, which is sensi- 
tive to disorder on the scale of the coherence length ~ 
50 A). We therefore expect the results to be representa- 
tive of the bulk. 

There have been several recent reports of measure- 
ments of A(r) for MgB2 by ac-susceptibility muon 
spin rotation and optical conductivity [Q techniques. 
These authors conclude that A(r) follows a power law 
dependence [A(r) - [||] and A(T) - T^ '^ g], which 
is at odds with tunneling and other measurements which 
indicate that there is a sizeable s-wave gap. Here we 
present the first high resolution measurements of A(r) in 
both single crystal and polycrystalline samples of MgB2. 
We find strong evidence for a predominately exponential 
temperature dependence of A at low temperature consis- 
tent with s-wave behavior. The gap deduced from fits to 
the data however, is significantly smaller than the BCS 
weak coupling value. 

Single crystal samples of MgB2 were prepared by a 
high pressure synthesis route as described in Ref. 
The data presented here was taken on a plate-like crys- 
tal with dimensions 0.35 x 0.22 x 0.1 mm'^, the smallest 
dimension being the c-axis. A second sample showed es- 
sentially identical behavior. The crystal orientation was 
verified with an x-ray Laue camera to within ^ 3°. The 
as-grown samples are shiny and gold in color but tarnish 
quickly in air. To ensure we measured clean surfaces we 



immersed the crystals briefly in a 0.5 % solution of HCl 
in cthanol, which removed this surface discoloration. 

Measurements were also made on polycrystalline MgB2 
samples, which were fabricated from commercially avail- 
able (Alfa Aesar) powder. This powder has a relatively 
wide range of grain sizes which complicates quantitative 
analysis of the data. For the analysis to give a reliable 
estimate of the absolute values of A(r) grain sizes ^ 5 A 
are required . To obtain such a distribution the powder 
was ground in an agate mortar and then sedimented in 
acetone for one hour The resulting powder was then 
cast in epoxy (~6% by volume). The grain size distri- 
bution was measured from scanning electron microscope 
images. The diameter of 96% of the grains was less than 
1/im, and the mean diameter was 0.56/im. 

Measurements of penetration depth were performed in 
a tunnel diode oscillator operating at 11.9 MHz p, with 
frequency stability of a few parts in 10^° Hz" 2. This 
translates to a resolution in A of 10~^^m for our powder 
samples and 10"^"m for the single crystals. A particu- 
lar feature of our apparatus is the very low value of the 
ac-probe field which we estimate to be ~ IfiT. Ambient 
dc fields are shielded to a similar level with a mu-metal 
can. Changes in the oscillator frequency are directly pro- 
portional to the inductance of the probe coil and hence 
to the susceptibility of the sample. For single crystal 
samples, this frequency shift can be directly related to 
changes in the penetration depth using the known sample 
dimensions ||^. For polycrystalline samples the relation 
between the measured susceptibility (per unit volume of 
superconductor) x ^nd A is more complicated and de- 
pends on the size distribution of the grains. For well 
separated grains 
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where Ni is the measured number of grains of radius 
ri ^ which are assumed to be spherical ||l^. A(T) of 
our polycrystalline sample was determined from the mea- 
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sured x(2^) ^^'^ grain size distribution by solving Eq. (|l|) 
at each temperature point. 

The superconducting transition of the single crystal is 
shown in the inset to Fig. 0. The Tc was 38.0 K and 
the 10-90% transition width was 0.3 K showing that the 
crystal was of high quality. The main part of this figure 
shows the low temperature (T <12 K) behavior of A(T) 
for this crystal with the field applied along each of the 
two principle crystallographic directions. AA(r) denotes 
the change in A relative to that at our base temperature 
(r=1.35 K). The striking feature of these curves is the 
lack of temperature dependence of AA(T) below around 
4 K, which points to the existence of an energy gap in all 
directions in fc-space. 

For H\\c only in-plane currents are excited and the 
in-plane penetration depth (Aab) is measured directly, 
within a calibration factor which depends on the dimen- 
sions of the crystal. This calibration factor was esti- 
mated using the procedure outlined in Ref. Q. As the 
faces of the crystal were somewhat uneven we estimate 
that the absolute values of AA(T) are only accurate to 
about 20%. For H\\ab both in-plane and out-of-plane 
currents screen the field and so we measure a mixture 
of Aab and the out-of-plane penetration depth Ac- As 
A ^ la,lb,lc ila,h,lc are the dimensions of the crystal 
along the a, b, and c directions respectively), the effective 
volume of the crystal penetrated by the field is approx- 
imately (for H\\a) V = 2a{lalb^b + lalc^c), where a is 
the field enhancement due to the demagnetising effects 
(a = 1.2 for _ff ||a6). In the figure we show an effective A, 
AAe = AV/{2alJb) = AAb + Ic/k^Xc, which would be 
equal to AA;, if the crystal were thin. Again there is some 
uncertainty in the absolute values (~ 20%) because of the 
uneven surfaces. For this reason we do not attempt to 
determine AAc from these data. However, we note that 
AAe(T') has a very similar T dependence to AAah but is 




FIG. 1. Temperature dependence of the penetration depth 
for an MgB2 single crystal. Data for H parallel to c (left axis) 
and ab (right axis) are shown. The solid lines are fits to Eq. 
(0). The inset shows the susceptibility near Tc- 



approximately 1.4 times larger. With the aspect ratio of 
the crystal {Ic/h) equal to 2.2, we conclude that AAc is 
between 1.5 and 2.5 times larger than AAab up to 12 K. 

For our crystals Aa6(0) ~ 1100 A, ^afc(O) ~55 A [|, 
and the mean free path (determined from the resistivity 
at 40 K) is - 250 A ||,|l|. We therefore expect to be 
in the moderately clean, local limit. Hence, for s-wave 
pairing we expect, for T ^ Tc/S, that A(T) should follow 
the BCS behavior 



AA(r) ~ Ao 



— ex I- — 
2T V T 



(2) 



were Aq is the value of the energy gap at zero temper- 
ature and Aq = A(0) in the weak coupling case. In Fig. 
|l| we show a fit of Eq. (|^) to the crystal data in the two 
orientations. The fit to AAafc {H\\c) gives Aq = 29±2 K, 
and to AAe {H\\ab) gives Aq = 32 ± 2 K. We note the 
remarkable similarity between these two values (the error 
reflects the different values of Aq obtained as the upper 
limit of the fit was varied from 10 K to 15 K). 

In Fig. H we show A(r) calculated from the susceptibil- 
ity of the sedimented polycrystalline sample. The data 
are very similar to that of the single crystals. A fit of the 
data to Eq. (||) gives Aq — 30±2 K, which is in very good 
agreement with the values found for the single crystal. 
The susceptibility (x) over the full temperature range is 
shown in the inset. The lack of sharp change in x at 
Tc for this sample, is entirely due to the small size of the 
grains [compared to A(T)] obtained by the sedimentation 
procedure and does not indicate inhomogeneity. This is 
confirmed by the calculated temperature dependence of 
the superfluid density (see below) and the sharp transi- 
tion observed in an unsedimented sample (not shown). 

As mentioned above, several other authors phH] have 
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FIG. 2. Temperature dependence of the penetration depth 
for a sedimented polycrystalline sample of MgB2 . Fits to the 
BCS expression [Eq. m] (solid line) and a law (dashed 
line) are shown. The inset shows the susceptibility over the 
full temperature range. 
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claimed that A(r) in MgB2 does not follow a simple expo- 
nential T dependence but rather a power law dependence 
with an exponent close to 2. In Fig. ^ we show a fit to 
the polycrystalline data. Clearly, this fit is much worse 
than the BCS dependence. We have also tried to fit other 
forms such as AT + BT'^ or AT'^/{T + T*) but no sig- 
nificant improvement was found. The same conclusion is 
reached from fits to the single crystal data. The key rea- 
son for this difference is that previous studies have been 
limited to temperatures above ~ 4 K, whereas AA(r) 
only shows a clear signature of gapped behavior below 
this temperature. 

From the polycrystalline data we are able to deduce 
absolute values of A(r). We find that A(0) = 1600 ± 
200 A which is in good agreement with other studies 
1^,^. In Fig. I^a we use this value to calculate the su- 
pcrfluid density p = [A(0)/A(T)]^ for this sample. This 
uncertainty in A(0) does not make any significant differ- 
ence to the temperature dependence of p{T). 

MgB2 has an anisotropic structure, with Mg atoms 
sandwiched between planar, hexagonal boron rings. It 
is therefore expected that there will be some anisotropy 
between the afo-plane and the c-axis responses. Recent 
results on aligned crystallites and single crystals have 
shown a significant anisotropy in Hc2 which implies an 
anisotropy in the coherence lengths, 7 — S,ab/S.c- In an 
anisotropic Ginzburg-Landau theory this implies a simi- 
lar anisotropy in A, 7 = Xc/Xat- There is some disagree- 
ment about the magnitude of this anisotropy, Ref. fl^ ] 
gives 7 = 1.7 ± 0.1 and Refs. ||,| give 7 = 2.6 ± 0.1. 
We note that the latter of these measurements were con- 
ducted on crystals identical to the ones measured here. 
As far as we know there is no general solution for the 
moment of a sphere when A is anisotropic, however so- 
lutions do exist in the limits X ^ r ||l^ and A <C r 
| p6[ . Fortunately, these two limits give similar results for 
7 < 3. To within ±10% the effective A(0) equals 1.2 A^t 
or 1.5 Xab for 7=1.6 and 2.6 respectively. We can use 
this fact to estimate the values of Xab{0) and Ac(0) from 
our polycrystalline data. We find that Aah=1300 A and 
Ac=2100 A or Xab ^ 1100 A and A^ ~ 2800 A , for 
the two values of 7 respectively. A value of Ac roughly 
two times Xab is consistent with the much stronger T de- 
pendence of AA in the H\\ab configuration as discussed 
above. 

Using the value of Xab = 1100 A (appropriate to the 
anisotropy of our crystals), we calculated the in-plane 
superfiuid density (Fig. ||b) from the single crystal data 
for AXab{T). The overall dependence is similar to that 
of the polycrystalline data, suggesting that there is not 
a significant anisotropy in the temperature dependence 
of p. Although the H\\ab data (AAg) is a mixture of Ac 
and Xab, we note that if we choose Ae(0) = 1750 A the 
calculated superfiuid density (pe = Ae(0)/Ae(T)) almost 
exactly overlaps Pab{T), showing that the two differ by 
only a simple scale factor and supporting the conclusion 



that there is little anisotropy in the temperature depen- 
dence of p. 

Although the overall temperature dependence of A de- 
duced from our measurements is consistent with an s- 
wave BCS picture, the value of the gap obtained is con- 
siderably smaller than the usual BCS weak coupling value 
(Ao/Tc = 1.76). Similar conclusions have been reached 
from some tunneling data fl^ , and specific heat studies 

The calculated band structure of MgB2 is composed 
of two sets of distinct bands; quasi- two-dimensional 
tubes which run along the c axis and a roughly three- 
dimensional network |l9[ . Several authors have proposed 
that there may be different gaps associated with these 
two bands [^,^. For these band specific gaps to re- 
main distinct the scattering between bands needs to be 
weak, and there needs to be some weak coupling between 
the bands so that the TcS are the same. Bouquet et al. 
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FIG. 3. (a) Superfiuid density [p = \\<d)/\\T)] of the 
polycrystalline sample (o), along with a fit to the multigap 
model (solid line). The contribution of the small gap (pA') 
and the large gap (p^b) in the model are also shown, (b) 
In-plane superfiuid densities of the single crystal sample [pab 
(o) and Pe (a)], along with a fit to the multigap model (solid 
line). The predicted behavior of pab and pc in the anistropic 
gap model for a = 2.2 is shown by the dashed lines (denoted 
pt^ and p^^). 
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||18| have recently proposed a simple phenomenological 
model to calculate the thermodynamic properties in such 
a two band system. The model is based on the so-called 
a model which is often used to describe strong coupling 
effects ||2l|] . It is supposed that the gaps on the two bands 
A'' and A'* follow the usual BCS T dependence, but have 
T = values which differ from the BCS one. The super- 
fluid densities from each band add together to give the 
measured total. These two gaps and the ratio of the 
contributions of each to the total superfluid density are 
fitting parameters. 

In Fig. ^ we show a least squares fit of both the poly- 
crystalline and single crystal data to this model. The 
agreement with the data is extremely good in both cases. 
We find that for the polycrystalline data, A* = 30 K, 
A^ = 89 K with relative proportion 40:60 respectively. 
For the in-plane single crystal data the parameters are 
found to be: A** = 29 K, A'' = 75 K with relative pro- 
portion 45:55. The small gap (A'*) is identical to that 
found from the fits of AA(r) to Eq. (|2|). Bouquet et al. 
have also applied this model to heat capacity data and 
find very similar results for both the gap values and their 
relative weights. 

First principles, density functional theory calculations 
of the electron phonon coupling in MgB2, have been per- 
formed by Liu et al. [Q. These calculations predict the 
existence of two distinct gaps with almost equal weights 
; a small gap (A/Tc — 0.65) associated with the 3D 
sheets and a larger gap {^/Tc = 2.0) associated with 
the 2D tubes. Both the size of the gaps and their rel- 
ative weights are in good agreement with the Bouquet 
model. In addition, because the smaller gap is located on 
the 3D sheet, we expect the effective gap, derived from 
fits to AA(r) at low temperature to be fairly isotropic. 
This is also in agreement with our observations. We note 
however, that Raman scattering measurements | |2^ on 
crystals prepared identically to those studied here, have 
clearly identified only one sharp pair-breaking peak at 
105 cm~^ corresponding to a gap value A(0)/rc = 2.0. 
Although there is some scattering below this peak, a sec- 
ond peak corresponding to a second smaller gap has not 
been seen. 

An alternative model |Q assumes that there is a single 
s-wave gap which is anisotropic, having a k dependence, 
A(z) = A(l -I- az'^)/{l + a), where z is the cosine of the 
polar angle and a is a fitting parameter which controls the 
anisotropy. Following the procedure outlined in Ref. ||2^ ] 
we have calculated the superfluid density in this model 
for arbitrary a. We find that to fit the observed ratio 
of A/Tc = 0.75 we need a = 2.2 ± 0.4, or a total gap 
anisotropy (1 + a) = 3.2 ±0.4. The overall T-dependence 
of pab calculated for this value of anisotropy is shown in 
Fig. |3|d and can be seen to be in serious disagreement 
with the data. Better agreement can be found if we set 
Aab(O) 650 A. However, this is almost a factor 2 below 
our estimate and well outside our expected error. The 



very different behavior of pc in this model (see figure) is 
also at odds with observed correspondence between pat 
and pe- Our data therefore clearly favor the multigap 
model above. 

In conclusion we have presented the first study of the 
penetration depth of both single crystal and polycrys- 
talline MgB2 samples. In agreement with other direct 
probes of the symmetry of the order parameter we find 
that A(r) is well described by an s-wave behavior, but 
with a minimum gap that is significantly less than the 
weak coupling BCS value. Of the various explanations 
for the small gap value we find that the two-gap model 
of Bouquet et al. is in best agreement with the data. We 
conclude therefore, that although the pairing interaction 
in this compound may be phononic, the gap structure is 
far from conventional. 
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F. Bouquet, K. Maki, P. Timms, J.F. Annett and J.R. 
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ported by the New Energy and Industrial Technology 
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